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/N i Following the prediction by Akaishi and Yamazaki of relatively narrow if-nuclear states, 

?*Zl ' deeply bound by over 100 MeV where the main decay channel KN — > 7rS is closed, sev- 

25T ' eral experimental signals in stopped K~ reactions on light nuclei have been interpreted 

recently as due to such states. In this talk I review (i) the evidence from K~ -atom data 
for a deep ^"-nucleus potential, as attractive as V^-(po) ~ —(150 — 200) MeV at nuclear 
matter density, that could support such states; and (ii) the theoretical arguments for a 
shallow potential, V R (po) ~ —(40 — 60) MeV. I then review a recent work by Mares, 
O , Friedman and Gal in which i?-nuclear bound states are generated dynamically across 

the periodic table, using a RMF Lagrangian that couples the K to the scalar and vector 
meson fields mediating the nuclear interactions. The reduced phase space available for 
K absorption from these bound states is taken into account by adding a density- and 
energy-dependent imaginary term, underlying the corresponding i?-nuclear level widths, 
K^ , with a strength constrained by K~ -atom fits. Substantial polarization of the core nu- 

» . ■ cleus is found for light nuclei, with central nuclear densities enhanced by almost a factor 

C^ ' of two. The binding energies and widths calculated in this dynamical model differ ap- 

preciably from those calculated for a static nucleus. These calculations provide a lower 
limit of Y R ~ 50 ± 10 MeV on the ' 
in the range B R = 100 - 200 MeV. 



1. Introduction 

The -ftT-nucleus interaction near threshold is strongly attractive and absorptive as 
suggested by fits to the strong- interaction shifts and widths in if^-atom levels 1 . 

Global fits yield 'deep' optical potentials V R (p ) (150 - 200) MeV 2 < 3 ' 4 * 5 , 

whereas other, more theoretically inclined works that fit the low-energy K~p reac- 
tion data, including the 1 = unstable bound state A(1405) as input for construct- 
ing a density dependent optical potential, suggest relatively 'shallow' potentials 

with V R (po) (40 - 60) MeV 6 > 7 > 8 > 9 . The issue of the depth of the attractive 

^-nucleus potential is briefly reviewd in Sections 2,3. 

* Dedicated to Walter Greiner on the occasion of his 70th birthday, presented at the International 
Symposium on Heavy Ion Physics, ISHIP 2006, Frankfurt, April 2006. 
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Fig. 1. Left: calculated spectrum of K 'deeply bound' atomic states in 208 Pb. Right: saturation 
of width r for the 2p K~- 208 Ph atomic state as function of Im bo, for Re bo = 0.62 fm. 



It is instructive at this introductory stage to mention that irrespective of the 
depth of the real part of the potential, and quite paradoxically due to its strong 
imaginary (absorptive) part, relatively narrow K~ deeply bound atomic states were 
predicted to exist 10 ' 11 , also confirmed in Ref. 8 . Figure 1 from Ref. n shows on 
the left-hand side a calculated spectrum of K~ atomic states in 208 Pb where, in 
particular, all the circular states below the 1% (I = 6) state are not populated by 
X-ray transitions due to the strong if _ -nuclear absorption, and on the right-hand 
side it demonstartes saturation of the 2p atomic-state width as a function of the 
absorptivity parameter Im &o of the optical potential. The physics behind is that a 
strong imaginary part of V^ acts repulsively, suppressing the atomic wavefunction in 
the region of overlap with Im V^. The calculated width of the 'deeply bound' atomic 
Is and 2p is less than 2 MeV, calling for experimental ingenuity how to form these 
levels selectively by a non radiative process 12 . Relatively narrow p deeply bound 
atomic states were also predicted, due to the same width saturation mechanism 11 . 

This saturation mechanism does not hold for nuclear states which retain very 
good overlap with the potential. Hence, the question to ask is whether it is possible 
at all to bind strongly K mesons in nuclei and are such bound states sufficiently 
narrow to allow observation and identification? The current experimental and the- 
oretical interest in this question was triggered back in 1999 by the suggestion of 
Kishimoto 13 to look for such states in the nuclear reaction (K~ ,p) in flight, and by 
Akaishi and Yamazaki 14 who suggested to look for a KNNN 1 = state bound 
by over 100 MeV for which the main KN — > 7r£ decay channel would be kincmati- 
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cally closed. Some experimental evidence has been presented recently for candidate 
states in (K~ top7 n) and (K^. op ,p) reactions on 4 He (KEK-PS E471, 15 ' 16 respec- 
tively), in the (K~ ,n) in-flight reaction on 16 (BNL-AGS, parasite E930 17 ), and 
in K~ t reactions in Li and 12 C, observing back-to-back Kp pairs from K~pp — ► Kp 
(DA$NE-FINUDA 18 ). None of it is sufficiently conclusive. 

It is interesting then to study theoretically K nuclear states in the range of bind- 
ing energy B^ ~ 100 — 200 MeV, as suggested by these recent experiments, and in 
particular the width anticipated for such deeply bound states. The relatively shal- 
low AT-nucleus potentials 8 ' 9 derived from the microscopic construction by Ramos 
and Oset 7 , which are discussed in Section 2, are of no use in this context, since 
they cannot generate even within a dynamical calculation binding energies sub- 
stantially greater than the potential depth of about 50 MeV. One must therefore 
depart from the microscopic approach in favor of a more phcnomcnologically in- 
clined model which is constrained by data other than two-body KN observables. 
The theoretical framework described here in Section 4, due to Mares, Friedman and 
Gal 5 ' 19 , is the relativistic mean field (RMF) model for a system of nucleons and one 
K meson interacting through the exchange of scalar (a) and vector (uJ, p) boson 
fields which are treated in the mean-field approximation. The RMF is a systematic 
approach used across the periodic table beyond the very light elements explored by 
other techniques, and it can be used also to study multi-A' configurations and to 
explore the K condensation limit 20 ' 21 . Similar RMF calculations have been recently 
reported by the Frankfurt group for A^ states in nuclei 22 ' 23 . 



2. -RT-nucleus potential from chirally motivated models 

The Born approximation for the if-nucleus potential due to the leading-order 
Tomozawa- Weinberg (TW) vector term of the chiral effective Lagrangian 24 , 

where /„■ ~ 93 MeV is the pseudoscalar meson decay constant, yields sizable at- 
traction Vx(po) ~ —55 MeV for po = 0.16 fm~ 3 . Iterating the TW term plus 
next-to-leading-order terms, within an in-medium coupled-channel approach con- 
strained by the KN — 7r£ — it A data near the KN threshold, roughly doubles this 
A'-nucleus attraction. It is found in these calculations (e.g. Ref. 25 ) that the A(1405) 
quickly dissolves in the nuclear medium at finite densities, well below po, so that 
the repulsive free-space scattering length a K - p becomes attractive, and together 
with the weakly density dependent attractive a K - n it yields an attractive density 
dependent effective isoscalar scattering length a c ff(p), denoted here also as &o(p) 
(positive for attraction): 



b a(p) = ^K-p(p) + a K - n {p)) , b (po) ~ 0.9 fm , (2) 
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Fig. 2. Re a B ff (top) and Im a a g (bottom) for the KN scattering length as function of the density, 
calculated without requiring self consistency (dashed) and requiring it (dot-dash and solid). 



leading to a strongly attractive if-nucleus optical potential: 

2tt 



VR(r) = - 



f-KN 



Mp)p{t) , Re V R {p ) ~ -110 MeV. 



(3) 



Here jjlkn is the KN reduced mass. However, when V^ is calculated self consistently 
(SC), namely by including V R in the in- medium propagator used in the Lippmann- 
Schwinger equation determining Vj( as a ^T-nucleus T matrix, the resultant K- 
nucleus potential is only moderatelly attractive, with depth between 40 — 60 MeV. 
This is shown in Fig. 2, where the reduction of Re a e ff from the (non SC) dashed 
line to the (SC) solid line is clearly seen 9 . The reason for this weakening of the 
T matrix, approximately back to its TW potential kernel Eq. (1), is the strong 
absorptive effect of V^ in the if-nuclcus propagator which suppresses the higher- 
order TW potential. It is the main reason that different theoretical calculations 
that obtain ^"-nucleus potential depths ~ 100 MeV without requiring SC, find 
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ultimately upon requiring SC depths ^ 50 ± 10 MeV at po 6 ' 7,9 - The Akaishi and 
Yamazaki scheme of calculation 14 differs in many respects from the theoretical 
derivations here outlined. It gets rid of the 7rS coupled-channel effects too early in 
the calculation, essentially using a one-channel KN energy-independent complex 
potential that fits the / = 0, 1 scattering lengths and the position of the A(1405) 
unstable bound state for I = 0. This kind of construction is far from being unique, 
certainly with respect to models that use a comprehensive set of KN low-energy 
data 7,9 . Furthermore, their analog of the TW potential kernel is considerably more 
attractive; and moreover, they do not impose the SC requirement, and it appears 
that the absorptive effect mentioned above is missing in their scheme. Therefore, 
their calculations have little common with the other 'theoretical' calculations in 
this field. 
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Fig. 3. Real part (left) and imaginary part (right) of the ftT-nucleus potential for 58 Ni and 208 Pb, 
obtained in a global fit to ftT~-atom data, for a tp potential, for the DD potential 1 and for potential 



see text. 



3. /^-nucleus potential from fits to K -atom data 

The X~-atom data used in global fits 1 span a range of nuclei from 7 Li to 238 U, with 
65 level-shift and -width data points. Figure 3 shows the real part and the imaginary 
part of the A'-nucleus potential for 58 Ni and 208 Pb for a tp potential, where the 
complex strength t is fitted to these data, and for two density-dependent potentials 
marked by DD and F, also fitted to the same data. For the real part, the depth of the 
tp potential for a typical medium- weight to heavy nucleus is about 80 MeV, whereas 
the density-dependent potentials are considerably deeper, 150 — 200 MeV. These 
latter potentials also yield substantially lower x 2 values of 103 and 84, respectively, 
than the value 129 for the tp potential. In particular, the shape of potential F 
departs appreciably from p(r) for p(r)/po < 0.2, where the physics of the A(1405) 
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still plays a major role. We note that, although the two density-dependent potentials 
have very different parameterizations, the resulting potentials are quite similar. In 
model F, considered in Ref. 5 , one divides configuration space into an 'internal' 
region and an 'external' region for both Re V K and Im V K : 

b a - B F(r) + b [1 - F{r)} , (4) 

using the function 

F(r) = — i— , x = (r- R x )/a x , R x = R^A 1 / 3 + 5 X . (5) 

e x + 1 

Here the parameter bo represents the interaction in the external region, and may be 
held fixed at its free KN value, whereas the parameter Bq represents the interaction 
inside the nucleus. The K~ atomic fit F used a fixed value a x = 0.4 fm, with 
little sensitivity for the precise value in the range 0.2 — 0.5 fm, optimizing for 
R x o = 1.30 ± 0.05 fm, 5 X = 0.8 ± 0.3 fm, thus implying that the modification of 
the free KN interaction, namely the transition from the bo term to the Bo term in 
Eq. (4), takes place at radii somewhat outside of the nuclear 'half-density' radius. 
Consequently, V K departs in shape from p{r) for densities less than about 20% of 
the central density. 

4. K"-nucleus RMF Dynamical Calculations 

4.1. K -nucleus RMF methodology 

In this model, expanded in Ref. 5 , the (anti)kaon interaction with the nuclear 
medium is incorporated by adding to Cn the Lagrangian density Ck 20,21 : 

C K = V'^KV^K - m K KK - g aK m K oKK . (6) 

The covariant derivative T>^ = d^ +ig UJ K l JJ^, describes the coupling of the (anti)kaon 
to the vector meson u. The coupling of the (anti)kaon to the isovector p meson was 
neglected, a good approximation for the light N = Z nuclei. Whereas extending the 
nuclear Lagrangian Cm by the Lagrangian Ck does not affect the original form of 
the corresponding Dirac equation for nucleons, the presence of K leads to additional 
source terms in the equations of motion for the meson fields a and loq to which the 
K couples. The K thus affects the scalar and vector potentials which enter the 
Dirac equation for nucleons, and this leads to a rearrangement or polarization of 
the nuclear core. 

In order to preserve the connection to studies of K~ atoms, the Klein Gordon 
(KG) equation of motion for the K is written in the form 9 : 

[A-2p(B s -v-+V K + V c ) + (V c + B s -i>-) 2 ]K = (h = c=l), (7) 

where V c is the K~ static Coulomb potential and /i is the if-nucleus reduced mass. 
The superscript s.p. in B sp - — B s ^'+iT K /2 stands for the single-particle K binding 
energy which is equal to the K separation energy only in the static calculation. The 
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difference B^' — B^ provides a measure of the nuclear rearrangement energy. The 
real part of the K optical potential V^ is then given by 

/i 2 rriK ZrriK 

where S = g a KO~ and V = g^K^o & re scalar and vector potentials due to the a and lo 
mean fields, respectively. Whereas Re Vjj is implicitly state dependent through the 
dynamical density dependence of the mean-field potentials S and V, as expected 
from a RMF calculation, it is here also explicitly state dependent through the 
[1 — {B S f^' + V r c )/?rix] energy-dependent factor multiplying the vector potential V. 

Since the RMF approach does not address the imaginary part of the potential, 
Im V op t was taken in a phenomenological tp form, where its depth was fitted to 
the K~ atomic data 4 with Im bo = 0.62 fm. Note that the density p in the present 
calculations is no longer a static nuclear density, but is a dynamical entity affected by 
the K interacting with the nucleons via boson fields. Suppression factors multiplying 
Im V ov t were introduced from phase-space considerations, taking into account the 
binding energy of the antikaon for the initial decaying state, and assuming two-body 
final-state kinematics for the decay products in the KN —* irY mesonic modes as 
well as in the KN — ► YN nonmesonic modes. 

The coupled system of equations for nucleons and for the electromagnetic vector 
field Aq, for the p meson mean field, and for the mean fields a and luq with terms 
explicitly due to K coupling, as well as the KG equation (7) for K~, were solved 
self-consistently using an iterative procedure. In order to produce different values 
of binding energies, g a K was scaled down successively from its initial value g^ m 
appropriate to the _ftT~-atom fit and, once it reached zero, g^K too was scaled down 
until the K~ Is state became unbound. Obviously, the good global fit to the atomic 
data is thus lost. Furthermore, in order to scan the region of large values of B K - , of 
order 200 MeV, we also scaled up g a K from its initial value, while keeping g U K fixed 
at g^K m - To give rough idea, whereas the static calculation gave B K - = 132 MeV 
for the K~ Is state in 12 C, using g^K m ^9aK m ^ the dynamical calculation gave 
B K - = 172 MeV for this same state. 



4.2. Results for binding energies and widths 

A comparison between the statically calculated (dashed lines) and the dynamically 
calculated (solid lines) B K - and T K - for the Is state in K i 2 C is shown in Fig. 4 as 
a function of the coupling-constant strenghts . We note that the width calculated 
dynamically for the Is nuclear state in ^-G does not fall below 50 MeV for the 
range of variation shown, whereas the corresponding limiting value of the statically 
calculated width is about 35 MeV. Another feature shown in Fig. 4 concerns the 
effect of the imaginary potential on the binding energy: the dynamically calculated 
binding energy Br- when Im V^ is switched off is shown by the dotted line. It is 
clear that the absorptive potential Im V^ acts repulsively and its inclusion leads to 
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Fig. 4. Is K binding energy and width in ^ C calculated statically (dashed lines) and dynam- 
ically (solid lines) for the nonlinear RMF model NL-SH 26 as function of a^ = g^nlO^K^ an< ^ 
a a = g&K/g^x™'- a Lj is varied in the left panels as indicated, with o CT = 0, and a a is varied in 
the right panels as indicated, with a u = 1. The dotted line shows the calculated binding energy 
when the absorptive K~ potential is switched off in the dynamical calculation. 



less binding, particularly at low binding energies. The repulsive effect of Im V^ gets 
weaker with B K - , along with the action of the kincmatical phase-space suppression 
factors, and beginning with B K - ~ 100 MeV it hardly matters for the calculation 
of B K - whether or not Im V/j is included. The A dependence of these features is 



considered in Rcf. 

with 12 C: 



leading to the following conclusions beginning approximately 



The K binding energy saturates, except for a small increase due to the 
Coulomb energy. 

The difference between the binding energies calculated dynamically and 
statically, B^ n — B s ^ at , is substantial in light nuclei, increasing with B^ 
for a given value of A, and decreasing monotonically with A for a given 
value of Bj^. It may be neglected only for very heavy nuclei. The same 
holds for the nuclear rearrangement energy B S ^P — B^ which is a fraction 



of B 



dyn 
K 



r?stat 

B R . 



The width T(B^) decreases monotonically with A, as shown in Fig. 5. 
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Fig. 5. Dynamically calculated widths of the Is K -nuclear state in „^?C, t?Ca and ^? 8 Pb 
as function of the K~ binding energy for nonlinear RMF models. The dotted line is for a static 
nuclear-matter calculation with po = 0.16 fm . 



Figure 5 shows calculated widths T K - as function of the binding energy B K - 
for Is states in ^-i 2 C and x l°Ca, using the nonlinear NL-SH version 26 of the RMF 



model, and in K 20S Pb using the NL-TM1 version 27 . The dotted line shows the static 



208-1 
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'nuclear-matter' limit 
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(9) 



where / is the phase-space suppression factor used in the dynamical calculations, 
and T K L is given by 



K- 



Hkn 



Im b p 



(10) 



for the static value Im bg 
Eq. (10) holds for a K~ 



0.62 fm used in the calculations and for po — 0.16 fm -3 . 
Schrocdinger wavefunction which is completely localized 
within the nuclear central-density po region. The additional factor (1 — B K - /tuk) 
in Eq. (9) follows from using the KG equation rather than the Schroedingcr 
equation. It is clearly seen that the dependence of the width of the K~ nuclear 
state on its binding energy follows the shape of the dotted line for the static 
nuclear-matter limit of T^--, Eq. (9). This dependence is due primarily to the 
binding-energy dependence of the suppression factor / which falls off rapidly un- 
til Bk- ~ 100 MeV, where the dominant KN — ► 7rE gets switched off, and then 
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stays rather flat in the range B K - ~ 100 — 200 MeV where the width is dom- 
inated by the two-nucleon absorption modes. The larger values of width for the 
lighter nuclei are due to the dynamical nature of the RMF calculation, whereby 
the nuclear density is increased by the polarization effect of the K~ . The widths 
calculated in the range Br-- ~ 100 — 200 MeV assume considerably larger values 
than what the atatic calculation of the dotted line shows (except for Pb in the 
range B K - ~ 100— 150 MeV). Adding perturbatively the residual width neglected 
in this calculation, mainly due to the KN — ► 7rA secondary mesonic decay channel, 
a representative value for a lower limit T^ ~ 50 ± 10 MeV holds in the binding 
energy range B K - ~ 100 - 200 MeV. 

4.3. Nuclear polarization effects 



0.20 r 




100 150 200 250 

B K - (MeV) 



Fig. 6. Calculated average nuclear density p for _ C, _ O, _ Ca and _ Pb as function of 
the Is K~ binding energy for the NL-SH RMF model 26 . 



Figure 6 shows the calculated average nuclear density p = -j J p 2 dv as a function 
of B K - for the same K~ nuclear Is states as in Fig. 5 and for Is states in ^-O. 
The average nuclear density p increases substantially in the light K~ nuclei, for the 
binding-energy range shown here, to values about 50% higher than for these nuclei 
in the absence of the K~ meson. The increase of the central nuclear densities is 
even bigger, as demonstrated in the next figure, but is confined to a small region 
of order 1 fm from the origin. 

Figure 7 shows on the left-hand side calculated nuclear densities of ^- 2 C for 
several values of the Is K~ binding energies. The purely nuclear density, in absence 
of the K~ meson, is given by the dashed curve. The maximal value of the nuclear 
density is increased by up to 75% in the range of binding energies spanned in the 
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Fig. 7. Calculated nuclear density p of „1?C (left) and of 3J. 8 Pb (right) for several B K - values of 
Is K~ nuclear states, using the nonlinear RMF models NL-SH 26 and NL-TM1 27 , respectively. 



figure, and the enhancement is close to uniform over the central region of r < 1 fm, 
decreasing gradually to zero by r — 2 fm which already marks the nuclear surface. 
In this fairly small nucleus, the density is enhanced over a substantial portion of the 
nucleus. This is different than in heavier nuclei, as shown on the right-hand side of 
Fig. 7 for ^°. 8 Pb. Here, the enhancement of the maximal density (at the center of 
the nucleus) is appreciable only for B K - values above 150 MeV, reaching a factor 
of two, but it subsides almost completely by r = 2 fm which is still well within the 
nuclear volume. As a result, the average nuclear density of ^_ Pb shown in Fig. 6 
is only weakly enhanced as function of Bx- ■ 

The upper and lower panels of Fig. 8 show the calculated nuclear rms radius 
and the Is and lp neutron single-particle energies E n , respectively, for K -0 as a 
function of B^- . It is clear that the polarization effect of the Is K~ bound state 
on the Is nuclear core is particularly strong. The differences between the linear and 
nonlinear models reflect the different nuclear compressibility and the somewhat 
different nuclear sizes obtained in the two models. The increase in the nuclear rms 
radius of K l -0 for large values of B K - is the result of the reduced binding energy 
of the lpi/2 state, due to the increased spin-orbit term. 



5. Conclusions 

In this talk I reviewed the phenomenological and theoretical evidence for a substan- 
tially attractive i^-nucleus interaction potential, from a 'shallow' potential of depth 
40 - 60 MeV to a 'deep' potential of depth 150 - 200 MeV at p Q . I then reported 
on recent dynamical calculations 5 of deeply bound K~ nuclear states across the 
periodic table, aimed at providing lower limits on the widths expected for binding 
energies in the range of 100 — 200 MeV. Substantial polarization of the core nucleus 
was found in light nuclei. Almost universal dependence of K widths on the binding 
energy was found, for a given nucleus, reflecting the phase-space suppression factor 
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Fig. 8. Nuclear rms radius and neutron single-particle energies for _ O as function of the Is K 
binding energy, for the linear RMF model L-HS 28 (open circles) and the nonlinear RMF model 
NL-SH 26 (solid circles). 



on top of the increase provided by the density of the compressed nuclear cores. 
The present results already provide useful guidance for the interpretation of recent 
experimental results by placing a lower limit T^ ~ 50 ±10 MeV. For lighter nuclear 
targets such as 4 He, where the RMF approach becomes unreliable but where nuclear 
polarization effects are found larger using few-body calculational methods 14 ' 29 , one 
anticipates larger widths for K deeply bound states, if such states do exist. 
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